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Abstract: A discussion of the time scales in the s process appears to be an approriate aspect to discuss 
at the occasion of Roberto's TO*** anniversary, the more as this subject has been repeatedly addressed 
during the 20 years of collaboration between Torino and Karlsruhe. The two chronometers presented in 
this text were selected to illustrate the intense mutual stimulation of both groups. Based on a reliable 
set of accurate stellar (n, 7) cross sections determined mostly at FZK, the Torino group succeeded 
to develop a comprehensive picture of the various s-process scenarios, which are most valuable for 
understanding the composition of the solar system as well as for the interpretation of an increasing 
number of astronomical observations. 
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1 Introduction 

The question of time scales is directly related to a 
number of key issues in s-process nucleosynthesis and 
from the very beginning the corresponding chronome- 
ters were considered to represent an important source 
of information. The long-lived radioactivities, which 
are produced in the s process, were studied first be- 
cause of their cosmological importance and since they 
are less sensitive to details of stellar scenarios. But in 
principle any unstable isotope along the s-process reac- 
tion path can be understood as a potential chronome- 
ter, provided that it defines a time scale corresponding 
to a significant quantity. 

The shortest possible time scale is related to con- 
vective mixing in He shell flashes, which take place in 
thermally pulsing low mass asymptotic giant branch 
(AGB) stars, where turnover times of the less than an 
hour are attained (Sec. 0). 

Neutron capture occurs on time scales of days to 
years, depending on the s-process scenario. The life 
time of a given isotope A is determined by the stellar 
neutron flux, n„ x vt, and by the stellar (n,7) cross 
section, 
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where n„ denotes the neutron density and vt the mean 
thermal velocity. If isotope A is unstable against (3- 
decay with a life time comparable to A„, the reac- 
tion path is split into a branching with a characteris- 
tic abundance pattern that reflects this time scale and 



provides a measure for the neutron density at the stel- 
lar site lfT3 : 143I). This situation can be complicated by 
the fact that in most cases only theoretical evaluations 
are available for the neutron capture cross sections of 
the unstable isotopes. Furthermore, the beta decay 
the /3-decay rate of the branch point isotope A may 
depend on temperature and/or electron density of the 
stellar plasma as in case of ^'^^Lu (Sec. O. 

If the neutron capture time is comparable to the 
duration of the neutron bursts, the abundance pattern 
of such a branching can be used to test the time scale 
of the He shell flashes in AGB stars. Suited branch- 
ings of this type could be those at *^Kr and ^^Zr with 
neutron capture times of a few years. In these cases, 
the neutron flux dies out before reaction equilibrium 
is achieved. However, both branchings are difBcult to 
analyze because of signiflcant contributions from the 
weak s-process component from massive stars or from 
the r process, and a constraining analysis is difficult 
to achieve. 

The transport to the stellar surface in the third 
dredge-up phase can be investigated by means of the 
observed Tc abundances (Mathews et al. 1986; Smith 
and Lambert 1988; Busso et al. 1995). Analyses 
of these observations have to consider that the ter- 
restrial decay rate of ®^Tc (ti/2 = 2.1 x 10^ years) 
is reduced to a few years at s-process temperatures 
US; 113)) which implies that it is quickly cooled to 
temperatures below 10* K after its production. A com- 
plementary, temperature-independent chronometer for 
the third dredge-up is ^''Zr, which can be followed by 
the appearance of its daughter ^^Nb ((25h . 
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The time scale for the formation of the solar sys- 
tem can, in principle, be inferred from the abundance 
patterns, which are affected by the decayof nuclei 
with half-lives between 10^ and 10^ years l|2ai . Quan- 
titative studies based on isotopic anomalies found in 
presolar grains have confirmed that such effects ex- 
ist for 26 Al, -^iCa, ^°Fe, ^^Zr, and ^"''Pd. A compre- 
hensive overview of this discussion was presented by 
Busso, Gallino, & Wasserburg (1999) (see also the con- 
tribution by R. Reifarth to this volume). Also ^''^Pb 
was discussed as a potentially promising chronometer 
(Yokoi, Takahashi, & Arnould 1985), whereas '^^Mn, 
^^''I, and ^^^Hf were found to result from the contin- 
uous pollution of the interstellar medium by explosive 
nucleosynthesis in supernovae (Busso et al. 1999). 

A number of attempts have been made to constrain 
the cosmic time scale by means of s-process abundance 
information. In the course of these studies it turned 
out that the half-life of the most promising case, ^^^Lu 
(Audouze, Fowler, & Schramm 1972; Arnould 1973), 
was strongly temperature-dependent, making it an s- 
process thermometer rather than a cosmic clock as dis- 
cussed in Sec. [S] The other long-lived species, ''"K ^ 
and *^Rb (Beer and Walter 1984) are produced by at 
least two different processes and are difficult to inter- 
pret quantitatively. Therefore, recent analyses of nu- 
clear chronometers for constraining the cosmic time 
scale concentrate on the r-process clocks related to 
the decay of the long-lived actinides (itI) and of ^^^Re 
(Yokoi et al. 1983; Arnould, Takahashi, & Yokoi 1984; 
Mosconi et al. 2007; see also the contribution by A. 
Mengoni to this volume). 

In the following sections we will focus on the short- 
est s-process time scale related to the fast convective 
mixing during the He shell fiashes in AGB stars. 



2 The branching at ^^^I 

Xenon is an element of considerable astrophysical in- 
terest. The origin of the lightest isotopes, '^^'^Xe and 
'^^^Xe, can be exclusively ascribed to the so-called p 
process in supernovae (Arnould and Goriely 2003) . Their 
relative isotopic abundances are important for testing 
p-process models describing the proton-rich side of the 
valley of stability. Concerning the s process, xenon 
belongs to the six elements with a pair of s-only iso- 
topes. In this case, the relevant nuclei are ^'^^Xe and 
'^^"Xe, both shielded against the decay chains from 
the r-process region by their stable Te isobars. The 
abundances of these isotopes define the strength of the 
branching in the s-process reaction chain illustrated in 
Fig. [T] Since the p-process components of these r- 
shielded nuclei do not exceed a few percent, they are 
commonly considered to be of pure s origin. On the 
neutron-rich side, ^^*Xe and ^'^'^Xe can be ascribed to 
the r process since the f3~ half life of ^'^''Xe is short 
enough to prevent any significant s-process contribu- 
tions. Hence, the Xe isotope chain carries signatures 
of all nucleosynthesis scenarios that contribute to the 
mass region of the heavy isotopes with A > 90, and of- 
fers the possibility to constrain the underlying models. 
A detailed description of the isotopic abundance pat- 



tern of xenon involves necessarily quantitative models 
for all these processes. 
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Figure 1: The s-process reaction path between Te 
and Xe. The isotopes ^^^Xe and ^^°Xe are shielded 
against r-process contributions by their stable Te 
isobars. In contrast to "'^'^''Xe, ^^^Xe is partly by- 
passed due to the branching at The branching 
at ^^^Te is negligible unless the temperature is low 
enough that ground state and isomer are not fully 
thermalized. The branching at ^^'^I is unique since 
it results from the competition between (3~ and 
electron capture decays and is, therefore, indepen- 
dent of the neutron flux. 



The combined strength of the branchings in the s- 
process chain at ^^''Te and ^^^I (Fig. [l| is defined by 
the relative abundances of the s-only isotopes '^^^Xe 
and '^^"Xe. Both branchings are expected to be com- 
parably weak since only a small part of the total re- 
action flow is bypassing ^'^^Xe. Therefore, the {a)Ns 
value for '^^^Xe, which is characteristic of the reaction 
flow, is slightly smaller than the one for i=^°Xe. Since 
the solar isotopic ratio of the s-only isotopes is well 
defined (Pepin, Becker, & Rider 1995), the {a)Ns dif- 
ference can be obtained by an accurate measurement 
of the cross section ratio (Reifarth et al. 2002) . 

While the first branching at ^^^Te is marginal be- 
cause the population of ground state and isomer is 
quickly thermalized in the hot stellar photon bath (Taka- 
hashi and Yokoi 1987), leading to a strong dominance 
of the /3-decay channel, the second branching at ^^^I 
is utterly interesting. In contrast to all other relevant 
cases, this branching is only defined by the competi- 
tion between j3~ and electron capture decay (Figure 
[!}. Both decay modes are sufficiently short-lived that 
the neutron capture channel and hence the influence 
of the stellar neutron flux is completely negligible. 

Since the electron capture rate is sensitive to tem- 
perature and electron density of the stellar plasma 
(40), this branching provides a unique possibility to 
constrain these parameters without interference from 
the neutron fiux. Under stellar conditions the corre- 
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spending branching ratio at I is 

/- = ^p-/i^l3- + ^pEc) = 1 — ^i3Ec/{X,3- + \pEc). 

While the /3~-rate varies only weakly, the electron cap- 
ture rate depends strongly on temperature due to the 
increasing degree of ionization. Furthermore, at high 
temperatures, when the ions are fully stripped, the 
EC rate becomes sensitive to the density in the stellar 
plasma via electron capture from the continuum (Table 
[!}. The relatively small change of the branching ratio 
did not permit quantitative analyses until the stellar 
(n,7) rates of the involved isotopes were measured to 
an accuracy of 1.5% (jsil). 

Table 1: Beta-decay branching ratio /_ at as 
a function of electron density n,, and temperature 

(BO). 
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2.1 Thermally pulsing AGB stars 

The s-process abundances in the mass range A > 90 
are produced during helium shell burning in thermally 
pulsing low mass AGB stars (|36h by the subsequent 
operation of two neutron sources. The ^^C(a,n)^^0 
reaction, which occurs under radiative conditions at 
low temperatures [kT ~8 keV) and neutron densities 
[un < 10^ cm~'^) between convective He-shell burning 
episodes, provides most of the neutron exposure. The 
resulting abundances are modified by a second burst 
of neutrons from the ^■^Ne(a, n)'^^Mg reaction, which is 
marginally activated during the highly convective He 
shell flashes, when peak neutron densities of n„ > 10^° 
cm~^ are reached at kT «23 keV. Although this sec- 
ond neutron burst accounts only for a few percent of 
the total neutron exposure, it is essential for adjusting 
the final abundance patterns of the s-process branch- 
ings. It is important to note that the (n, 7) cross sec- 
tions in Te-I-Xe region are large enough that typical 
neutron capture times are significantly shorter than 
the duration of the two neutron exposures. Therefore, 
the abundances can follow the time variations of the 
neutron density. 

2.2 Convection in He shell flashes 

The effect of convection for the branchings at ^ = 
127/128 was extensively studied by means of the stellar 
evolution code FRANEC (Straniero et al. 1997; Chi- 
efli and Straniero 1989) using a time-dependent mixing 
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Figure 2: The calculated convective velocities as a 
function of the internal radius for 5 models (t=0, 
0.31, 0.96, 2.42 and 5.72 yr after the maximum of 
the TP) in a 3 Mq star of solar composition. The 
convective turnover time is a few hours only. The 
scale on the abscissa starts at the bottom of the 
convective shell. 



algorithm to treat the short time scales during ther- 
mal pulses along the AGB. Convective velocities were 
evaluated by means of the mixing length theory where 
the mixing length parameter was calibrated by a fit of 
the solar radius. An example for these results is given 
in Fig. [21 which represents a typical thermal pulse for 
a 3 Mq AGB star of solar composition. The calculated 
convective velocities are plotted for five models (t=0, 
0.31, 0.96, 2.42 and 5.72 yr after the maximum of the 
TP) as a function of the internal radius to show how 
the convective shell expands after the pulse maximum. 

The rather short convective turnover times of less 
than one hour that can be derived from Fig. [5] are in 
any case, shorter than the time during which the tem- 
perature at the bottom of the convective pulse remains 
higher than 2.5 x 10* K. Hence, the crucial transport 
time from the hot synthesis zone to cooler layers is only 
of the order of minutes. This can have an impact on 
those potential s-process branchings, which are char- 
acterized by branch-point isotopes with strong thermal 
enhancements of their decay rates. Even if such half- 
lives are reduced to a few minutes at the bottom of 
the He shell flash, the signature of the branching can 
survive by the rapid mixing of processed material into 
the cooler outer layers of the convective zone. 

2.3 Branching analysis at A=127/128 

The (n, 7) cross sections required for describing the 
time evolution of the ^^®Xe and ^^"Xe abundances dur- 
ing the He shell flashes have been accurately measured 
(|3ll : [30) . This information is crucial in view of the fact 
that only a comparably small fraction of the reaction 
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flow is bypassing ^^*Xe (Table [T|. With the Xe cross 
sections compiled by Bao et al. (2000) for example, 
the ^■^^Xe/^'^'^Xe ratio could only be calculated with 
an uncertainty of about 20%, whereas this uncertainty 
was reduced to ±1.5% with the improved cross section 
data. Another source of uncertainty in the nuclear in- 
put is due to the theoretical /3-decay rates listed in 
Table [T] Variation of these rates by a factor of two 
affects the branching ratio by less than 2%, because 
the decay of is dominated by the /3~-mode. 

In fact, the branching analysis indicates a more 
complex situation than one might expect from the trends 
of the branching factor in Table[T] During the low tem- 
perature phase between He shell flashes, the neutron 
density produced via the ^^C(a, n)^''0 reaction is less 
than 10^ cm~'^. The branching at ^'^'^Te is completely 
closed, resulting in a ^^^Xe/'^^'^Xe abundance ratio of 
0.93 relative to the solar values at the end of the low 
temperature phase due to the effect of the ^'^^I branch- 
ing. After the onset of convection at the beginning of 
the He shell flash the production factors of ^^*Xe and 
^^°Xe differ by 8%, corresponding to the solar ratio, 
but are then modifled during the flash. Since the (n, 7) 
cross sections of both isotopes are large enough for 
achieving local {a)Ns equilibrium, one would expect 
that the production factor of ^'^^Xe quickly approaches 
that of ^=*°Xe. In this case, the final abundance ratio 
would clearly exceed the solar value. 

In contrast, one finds that the branching still exists 
even during the high temperature phase of the He shell 
flash. There are essentially two effects, which concur 
to explain this behavior: 

• During the peak of temperature and neutron 
density, the electron densities at the bottom of 
the convective He shell flash, i.e. in the s-process 
zone, are between 15 x lO^'^ cm"^ and 20 x 10^*^ 
cm~^. Therefore, the branching at ^^®I is never 
completely closed. Even at the peak tempera- 
tures of the He-shell flash, typically 3% of the 
flow are bypassing ^^*Xe. 

• Around the maximum of the neutron density, 
the branching at ^^'^Te is no longer negligible 
and leads to an additional decrease of the ^^*Xe 
abundance. 

In view of these effects the ^^^Te branching can be 
fully understood only if the stellar reaction network 
is calculated with sufficiently small time steps so that 
the time scale of convective turnover is properly con- 
sidered. If the neutron density is followed in time steps 
of 10^ s up to freeze-out one finds an enhancement of 
the '^^'^Te branching of 30% compared to calculations 
performed on a time grid of 10^ s. However, reducing 
the time steps to 3 x 10'' s had no further effect on the 
abundance ratios. 

If averaged over the AGB evolution, the stellar 
models yield a final abundance ratio of '^^Xe/^^^Xe 
= 0.466±0.015, consistent with the smaller value mea- 
sured in the solar system = 0.510±0.005 if an 8% p 
contribution to the solar ^'^^Xe abundance is taken into 
account. 

In summary, the abundance ratio of ^^*Xe and 
^^°Xe in the solar system could be eventually repro- 
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Figure 3: The s-process reaction flow in the Lu 
region. The strength of the lines indicates that 
neutron captures on ^^^Lu leading to the isomeric 
state in ^^^Lu are more probable than those to the 
ground state. 



duced by combining the effects of mass density, tem- 
perature, neutron density, and convective turnover in 
a consistent way. This success represents an impres- 
sive confirmation of the stellar s-process model related 
to thermally pulsing low mass AGB stars. 



3 The i^^Lu/i^ejjf pair 

The mass region of the rare earth elements (REE) rep- 
resents an important test ground for s-process models 
because the relative REE abundances are known to 
better than ±2%, which implies that the branchings 
in this region are reliably defined by the respective s- 
only isotopes. Systematic analyses of these branchings 
contribute essentially to the quantitative picture of the 
main s-process component that has been achieved with 
stellar models for thermally pulsing low mass AGB 
stars (GaUino et al. 1998; Busso et al. 1999). Also 
the stellar (n, 7) cross sections (except for the lightest 
REE) are large enough that the assumption of reac- 
tion flow equilibrium is well justifled with respect to 
the analysis of several prominent s-process branchings 
in this mass region. This means that the flnal abun- 
dances are essentially determined during the freeze-out 
phase at the decline of neutron density towards the end 
of the shell flash. 

Among the REE branchings, ^^''Lu is especially at- 
tractive for the intricate way in which nuclear physics 
can affect the actual s-process production yields. This 
is illustrated in Fig. |3] showing that the reaction path 
in the vicinity of lutetium is determined not only by 
the stellar neutron capture cross sections of "^'^'^''Lu 
and ^^*Hf, but also by the thermal coupling of isomer 
and ground state in ^^®Lu. 

Due to its long half-life of 37.5 Gyr, ^^'^Lu was ini- 
tially considered as a potential nuclear chronometer for 
the age of the s elements (Q). This possibility is based 
on the fact that ^^^Lu as well as its daughter ^^®Hf 
are of pure s-process origin, since both are shielded 
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against the r-process beta decay chains by their sta- 
ble isobar ^^'^Yb. In a straightforward approach, the 
reaction flow in the branching at y4 = 176 and, there- 
fore, the surviving s abundance of ^^^Lu and of '^^^Hf 
would be determined by the partial (n, 7) cross sec- 
tions of ^^^Lu feeding the ground and isomeric state 
of ^''^Lu. Since transitions between the two states are 
highly forbidden by selection rules, both states could 
be considered as separate species in the description of 
the s-process branching at ^4 = 176 (Fig. Q. While 
the isomer rapidly decays quickly (ti/2 ~ 3.7 h) by 
/^-transitions to ^''^Hf, the effective s-process yield of 
'^^^Lu appeared well defined by the partial cross sec- 
tion to ground state, thus providing an estimate the 
age of the s elements by comparison with the actual 
solar system value. 

However, Ward and Fowler (1980) noted that ground 
state and isomer of ^^^Lu are most likely connected by 
nuclear excitations in the hot stellar photon bath, since 
thermal photons at s-process temperatures are ener- 
getic enough to populate higher lying states, which 
can decay to the long-lived ground state and to the 
short-lived 123 keV isomer as well. In this way, the 
strict forbiddeness of direct transitions between both 
states is circumvented, dramatically reducing the ef- 
fective half-life to a few hours. As a result, most of 
the reaction flow could have been directly diverted to 
'^^^Hf, resulting in a ^^'^Lu abundance much smaller 
than observed in the solar system. That this temper- 
ature dependence had indeed affected the '^'^^Lu/^^^Hf 
abundance ratio was conflrmed soon thereafter by Beer 
et al. (1981, 1984), thus changing its interpretation 
from a potential chronometer into a sensitive s-process 
thermometer. 

The temperature dependence was eventually quan- 
tified on the basis of a comprehensive investigation of 
the level structure of ^'"'^Lu (Klay et al. 1991b; Klay 
et al. 1991a; DoU et al. 1999). The lowest mediating 
state was identified at an excitation energy of 838.6 
keV, which implies that thermally induced transitions 
become effective at temperatures above Ts = 1.5 — 2 
(where Tg is the temperature in units of 10* K). Ac- 
cordingly, ground state and isomer can be treated as 
separate species only at lower temperatures. In ther- 
mally pulsing low mass AGB stars this is the case be- 
tween convective He shell flashes when the neutron 
production is provided by the ^^C(a,n)^^0 reaction 
in the so-called ^^C pocket. Under these conditions 
the abundance of ^^''Lu is directly determined by the 
partial cross sections populating ground state and iso- 
mer. 

During the He shell flashes, the higher tempera- 
tures at the bottom of the convective region lead to 
the activation of the ^^Ne{a,n)^^Mg reaction. It is in 
this regime that the initial population of ground state 
and isomer starts to be changed by thermally induced 
transitions. This affects the ^'^^Lu^ in the long-lived 
ground state that is actually produced by neutron cap- 
ture in the burning zone as well as the Lu fraction cir- 
culating in the convective He shell flash. The latter 
part is exposed to the high bottom temperature only 
for rather short times and, therefore, less affected by 
the temperature dependence of the half life. Once pro- 



duced, by far most of the long-lived ^'^''Lu^ survives in 
the cooler layers outside of the actual burning zone. 

3.1 Nuclear physics input 

The stellar neutron capture rates for describing the 
reaction flow in Fig. [3] were determined by accurate 
time-of-fiight (TOF) measurements of the total (n, 7) 
cross sections for ^'^^Lu and ^^^Lu using a 47r BaF2 ar- 
ray (Wisshak et al. 2006). The Maxwellian averaged 
cross sections (MACS) deduced from these data are 
five times more accurate than the values listed in the 
compilation of Bao et al. (2000). However, neutron 
captures on ^^^Lu may feed either the ground state 
or the isomer in ^^^Lu. Therefore, the total (n, 7) 
cross section was complemented by a measurement of 
at least one of the two partial cross sections (cr|, a™). 
Since the corresponding reaction channels could not 
be distinguished in the TOF measurement Q), the 
activation technique was used to determine the partial 
cross section to the isomeric state in ^^''Lu at thermal 
energies of kT = 5.1 and 25 keV. 

Activation measurements of the partial (n, 7) cross 
section to the isomeric state in ^'^^Lu at or near a neu- 
tron energy of 25 keV were performed via the '^Li(p, n)'^Be 
reaction (Beer and Kappeler 1980; Allen, Lowenthal, 
& de Laeter 1981; Zhao and Kappeler 1991) and also 
using a filtered neutron beam from a nuclear reactor 
(Stecher-Rasmussen et al. 1988). In all experiments 
the induced activity after irradiation was detected via 
the 88 keV 7-transition in the decay of ^^^Lu"^. These 
data were recently complemented by an activation at 
kT = 5 keV, adapted to the lower temperature in the 
^^C pocket Iplb . where most of the neutron exposure 
is provided by the ^^C(a, n)^''0 source, which operates 
at fcT = 8 keV. 

The half- life of ^^'^Lu was shown to decrease drasti- 
cally at the temperatures reached during He shell burn- 
ing in thermally pulsing low mass AGB stars (Klay et 
al. 1991a; Doll et al. 1999). Contrary to the situ- 
ation at low excitation energy, where transitions be- 
tween ground state and isomer are strictly forbidden 
by selection rules, interactions with the hot stellar pho- 
ton bath lead to induced transitions to higher lying 
nuclear states, which can decay into the ground state 
and into the isomer as well (Fig. 

Up to Ts = 2 the mediating states can not be 
reached and the reaction path of Fig. [3] is completely 
defined by the partial capture cross sections feeding 
ground state and isomer in ^'^^Lu. This situation pre- 
vails in the ^'^C pocket between He shell flashes. Be- 
tween Tg = 2.2 and 3.0, ground state and isomer start 
to be increasingly coupled. At flrst, this coupling leads 
to an increasing population of the ground state, be- 
cause of the smaller energy difference between isomer 
and mediating state. It is due to this effect that more 
^'^^Lu is observed in nature than would be created in 
a "cool" environment with Tg < 2. In this regime, 
internal transitions, /3-decays, and neutron captures 
are equally important and have to be properly consid- 
ered during He shell flashes, where the final abundance 
pattern of the s-process branchings are established by 
the marginal activation of the ^^Ne(Q, n)^^Mg neutron 
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mediating state (ij) 




Figure 4: Schematic level scheme of ^^^Lu illus- 
trating the thermal coupling between the long- 
lived ground state and the short-lived isomer. 
While direct transitions between these states are 
strictly forbidden by selection rules, thermal ex- 
citations in the hot stellar photon bath populate 
a higher lying state with intermediate quantum 
numbers that can decay both ways. At first this 
link depopulates the isomer towards the ground 
state because it is more easily reached from the iso- 
mer. However, at sufficiently high temperatures, 
induced transitions from the ground state to the 
mediating state can also feed the isomer, result- 
ing in the destruction of ^"^^Lu via /3 decay of the 
short-lived isomer. The efficiency of this mecha- 
nism is determined by the lifetime of the medi- 
ating state and by the branching ratio V for decays 
towards the isomer. 

source reaction at thermal energies around kT = 23 
keV. 

3.2 Branching factor 

The branching factor /„ describing the split of the re- 
action flow at ^'^^Lu (Fig. [S]) can be expressed in terms 
of the neutron capture rate A„ — {fJi7e)vTnn (where 
CT176 denotes the MACS of ^'^''Lu), vt the mean ther- 
mal neutron velocity, and n„ the neutron density) and 
of the temperature and neutron density dependent /3- 
decay rate of ^'^''Lu. Because of the thermal efi^ects on 
the lifetime of ^'^^Lu sketched above, f„ = A„/(A„+A/3) 
becomes a complex function of temperature and neu- 
tron density (Klay et al. 1991a; DoU et al. 1999). 

The bottom temperature in the He flashes reaches 
a maximum at the maximum expansion of the con- 
vection zone and then declines rapidly as the convec- 
tive thermal instability shrinks. According to current 
AGB models of low mass stars, the neutron burst re- 
leased by the ^'^Ne neutron source lasts for about six 
years, as long as Tg > 2.5, while the convective in- 
stability lasts for a much longer period of about 150 
years. The maximum temperature in the He flash in- 
creases slightly with pulse number, reaching Tg ~ 3 in 



the more advanced pulses |3^. 

Although the second neutron burst during the He 
shell flashes contributes only a few percent to the total 
neutron exposure, the final composition of the Lu/Hf 
branching is, in fact, established during the freeze-out 
phase at the end of the thermal pulse. It is only dur- 
ing the high-temperature phase that the thermal cou- 
pling between ground state and isomer takes place. In 
this context, it is also important to remember that 
the (n,7) cross sections in the Lu/Hf region are large 
enough that typical neutron capture times are signifi- 
cantly shorter than the duration of the neutron expo- 
sure during the He shell flash. 

Previous s-process calculations for thermally puls- 
ing low mass AGB stars (Gallino et al. 1988; Gallino 
et al. 1998), which were performed by post-processing 
using full stellar evolutionary models obtained with the 
FRANEC code (Chiefli and Straniero 1989; Straniero 
et al. 1995; Straniero et al. 2003; Straniero, Gallino, 
& Cristallo 2006) were successful in describing the so- 
lar main s-process component fairly well (Gallino et 
al. 1998; Arlandini et al. 1999). For the treatment of 
the branching at ^'^^Lu, however, this approach had to 
be reflned by subdividing the convective region of the 
He shell flashes into 30 meshes in order to follow the 
neutron density proflle in sufficient detail and to take 
the strong temperature dependence of the branching 
factor /„ properly into account. The production and 
decay of ^'^^Lu was calculated in each mesh to obtain 
an accurate description of the flnal ^^^Lu/^'^^Hf ratio 
in the He shell flash. 

After each time step of less than one hour, which 
corresponds to the typical turnover time for convec- 
tion in the He shell flash (H^), the abundances from 
all zones were averaged in order to account for the 
fast convective mixing. This treatment is well justified 
because there is no efficient coupling between ground 
state and isomer at temperatures below Tg < 2.5, 
which holds for all meshes except the first few from the 
bottom. Interestingly, this is the temperature when 
neutron production via the ^■^Ne(a, n)'^^Mg reaction 
diminishes. 

During the main neutron exposure provided by the 
^^C(a, nY^O reaction thermal effects are negligible for 
the production of ^'^^Lu because the temperatures in 
the ^^C pocket of Tg « 1 are not sufficient to reach 
the mediating state at 838.6 keV. Accordingly, the 
^'^''Lu/^^^'Hf ratio is simply determined by the par- 
tial (n, 7) cross section to the ground state, which 
would always produce too much ^'^''Hf and too little 
^'^''Lu^. This situation is illustrated in Fig. [5] show- 
ing the relative production factors of '^^^Lu and ^^^Hf 
during the 15"" He shell flash of the AGB model with 
initial mass 1.5 Mq and half-solar metallicity, which is 
representative of the overall relative abundance distri- 
bution ejected during the AGB phase. The s-process 
abundances are given as number fractions Yi — Xi/Ai 
{X being the respective mass fractions), normalized to 
that of ^'^"Sm at the end of the He shell ffash. This 
isotope was chosen as a reference because it is the best 
known case of an unbranched s-only isotope among the 
REE. 

After the s-process nuclides synthesized in the ^^C 
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Figure 5: Evolution of the ^^^Lu (squares) and 
^^®Hf (circles) abundances during the 15*'' He shell 
flash in an AGB star with 1.5 Mq. The time scale 
starts when the temperature at the bottom of the 
convective shell reaches Ts — 2.5, i.e. at the onset 
of the ^^Ne(a, n)^^Mg reaction. The values are 
given as number fractions normalized to that of 
i^°Sm at the end of the He shell flash. 



pocket are engulfed by the convective He shell, one 
finds the ratio shown at t = in Fig. (5] which marks 
the onset of s process nucleosynthesis in the He shell 
flash. In this phase temperatures are high enough to 
populate the the mediating state at 838 keV, leading 
to a strong increase in the production of the long-lived 
ground state of ^''^Lu'^. As shown in Fig. [5] ^''^Hf 
is almost completely destroyed at the beginning of the 
flash, when neutron density and temperature are high- 
est and ^^^Hf is efficiently bypassed by the reaction 
flow. Correspondingly, the ^^®Lu abundance reaches a 
pronounced maximum. As temperature and neutron 
density decline with time, the branching towards ^^^Hf 
is more and more restored, but the final abundance 
at the end of the He shell flash remains significantly 
lower than at the beginning. Because the thermal cou- 
pling between ground state and isomer depends so crit- 
ically on temperature, the full s-process network was 
followed in each of the 30 meshes of the convective 



and temperature by the refined zoning of the convec- 
tive He shell flash. Only in this way overproduction 
factors could be obtained, which were compatible with 
the limits defined by the nuclear physics uncertain- 
ties and by the ^^^Lu decay. In contrast to a previous 
study (0) , the meanwhile improved cross section infor- 
mation together with the multi-zone approach for the 
He shell flash appears to settle the ^^^Lu puzzle. While 
the '^^^Hf abundance was before overproduced by 10 — 
20%, the present calculations yield final ^'^^Lu^ and 
^'^^Hf abundances (relative to solar system values) of 
1.04 and 0.96. These results do not yet consider the 
decay of the produced ^'^^Lu^, a correction that brings 
the two numbers even into closer agreement. 

Recently, additional information from Coulomb ex- 
citation measurements pdl) and photoactivation stud- 
ies |23 ) have triggered renewed interest in the Lu/Hf 
problem and seem to provide further constraints 
for the convection at the bottom of the He shell flash. 



4 Conclusions 

The elemental and isotopic abundances found in na- 
ture carry the signatures of their nucleosynthetic origin 
as well as of the related time scales. The branchings 
at ^^*I and ^'^^Lu discussed in Sees. [2] and [3] are char- 
acterized by the rather short time scale of convection 
in He shell flashes of low mass AGB stars. Examples 
dealing with much longer time scales are discussed by 
R. Reifarth and A. Mengoni in their contributions to 
this volume. 

To decipher the information contained in abun- 
dance patterns can be an intricate and complex pro- 
cess, which requires 

- accurate nuclear physics data, in particular the rates 
for production and transmutation by nuclear reactions, 
but often also nuclear structure properties for describ- 
ing the impact of temperature, and 

- sufficiently detailed astrophysical prescriptions for a 
realistic or at least plausible modeling of the relevant 
phenomena. It is this duality that we had the pleasure 
to pursue in the Torino-FZK liaison for more than two 
decades, and which, in return, was found to provide 
surprising constraints for the physics of stellar evolu- 
tion. 



Plausible solutions of the Lu/Hf puzzle are charac- 
terized by equal overproduction factors for '^'^^Lu and 
^'^^Hf, at least within the experimental uncertainties of 
the nuclear input data. The decay of long-lived ^^^Lu 
ground state in the interstellar medium prior to the 
formation of the solar system reduces the calculated 
overproduction factor of ^^^Lu by about 10%, leaving 
the more abundant ^^^Hf essentially unchanged. In 
the light of these considerations, overproduction fac- 
tors between 1.00 and 1.10 are acceptable for ^^^Lu 
and between 0.95 and 1.05 for ^''^Hf 

In the course of these investigations it turned out 
that the observed ^^^''Lu/^^^'Hf ratio could only be re- 
produced if the thermal coupling of ground state and 
isomer in ^^'^Lu in combination with the neutron den- 
sity was followed within the gradients of time, mass. 
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